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5^ | Some recently discovered nonperturbative strong-field effects in tensor-scalar 

theories of gravitation are interpreted as a scalar analog of ferromagnetism: 
"spontaneous scalarization" . This phenomenon leads to very significant de- 
^ . viations from general relativity in conditions involving strong gravitational 

fields, notably binary-pulsar experiments. Contrary to solar-system experi- 
ments, these deviations do not necessarily vanish when the weak-field scalar 
coupling tends to zero. We compute the scalar "form factors" measuring these 
deviations, and notably a parameter entering the pulsar timing observable 7 
through scalar-field- induced variations of the inertia moment of the pulsar. 
An exploratory investigation of the confrontation between tensor-scalar the- 
ories and binary-pulsar experiments shows that nonperturbative scalar field 
effects are already very tightly constrained by published data on three binary- 
pulsar systems. We contrast the probing power of pulsar experiments with 
that of solar-system ones by plotting the regions they exclude in a generic 
two-dimensional plane of tensor-scalar theories. 
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I. INTRODUCTION 



Einstein's general relativity theory postulates that gravity is mediated only by a long- 
range tensor field. It has been repeatedly pointed out over the years (starting with Kaluza 
U) that unified theories naturally give rise to long-range scalar fields coupled to matter with 
gravitational strength. This led many authors, notably Jordan ||, Fierz ||, and Brans and 
Dicke 0], to study, as most natural alternatives to general relativity, tensor-scalar theories 
in which gravity is mediated in part by a long-range scalar field. The motivation for such 
theories has been recently revived by string theory which contains massless scalars in its 
gravitational sector (notably the model- independent dilaton). 

We shall consider tensor-scalar gravitation theories containing only one scalar field, as- 
sumed to couple to the trace of the energy-momentum tensor. The simplest example of 
such a theory is a scalar field only coupled to the gravitational sector through a nonminimal 
coupling £i?$ 2 (see Section [VT] below). For a study of the observable consequences of general 
tensor-scalar theories (containing one or several scalar fields), see Ref. 0. 

Actually, one generically expects scalar fields not to couple exactly to the mass but to 
exhibit some "composition dependence" in their couplings to matter. However, a recent 
study of a large class of viable string-inspired tensor-scalar models |J has found that the 
composition-dependent effects represent only a very small fraction (~ 10~ 5 ) of the effec- 
tive coupling to matter. Such fractionally small composition-dependent effects would be 
negligible in the gravitational physics of neutron stars that we consider here. 

The most general theory describing a mass-coupled long-range scalar contains one arbi- 
trary "coupling function" A(ip) The action defining the theory reads 

+S m {i, m ;A'( v ) g ;j . (1.1) 

Here, G* denotes a bare gravitational coupling constant, i?* = gi^R*^ the curvature scalar 
of the "Einstein metric" g^ u describing the pure spin-2 excitations, and ip our long-range 
scalar field describing spin-0 excitations. [We use the signature -+++ and the notation 
g* = — det g^ u ] The last term in Eq. (|1.1|) denotes the action of matter, which is a functional 
of some matter variables (collectively denoted by ip m ) and of the "physical metric" g^ v = 
A 2 (ip)g^ v . Laboratory clocks and rods measure the metric g^ u which, in the model considered 
here, is universally coupled to matter. The reader will find in Eqs. (|6.1|)-(|6.7|) below an 



explicit example (nonminimally coupled scalar field) of how an action of the type ( |1.1D , 
involving two conformally related metrics g* u and g^ v = A 2 (ip)g^ u , can naturally arise. 

The field equations of the theory are most simply formulated in terms of the pure-spin 
variables (g^ u ,{p). Varying the action ( |1.1|) yields 

Svra / _ 1 
c 4 V 2 



R% = 2d^d v( p + T; v - -T* g ; u , (1.2a) 



u 9^ = -^«(^)^ , (1.2b) 



with T% v = 2c g* 5S m /5g^ u denoting the material stress-energy tensor in "Einstein units", 
and a(<f) the logarithmic derivative of A(ip) : 



2 



(1.3) 



[All tensorial operations in Eqs. ( |1.2| ) are performed by using the Einstein metric g*, e.g. 
= g"^V*V*, T* = g*T£ v .\ As is clear from Eq. ( |1.2b| ), the quantity a((p) plays the role 



of measuring the (field-dependent) coupling strength between the scalar field and matter. It 
has been shown in Refs. |||7j that all weak-field ( "post-Newtonian" ) deviations from general 
relativity (of any post-Newtonian order) can be expressed in terms of the asymptotic value 
of a (<p) at spatial infinity and of its successive scalar-field derivatives. Let <po denote the 
asymptotic value of ip at spatial infinity, i.e., the "vacuum expectation value" of ip far away 
from the considered gravitating system. Let us also denote: «o = a(<po), Po = da(ipo)/dipo, 
(3'q = dfl(<p )/d<p . At the first post-Newtonian approximation, deviations from general 
relativity are proportional to the Eddington parameters 

7 = 7Edd - 1 = -2ag/(l + al) , (1.4a) 

J = /3 Edd - 1 = i/3 ^/(l + a 2 ) 2 , (1.4b) 

while at the second post-Newtonian approximation there enters, beyond 7 and (3, two new 
parameters |||7| 

e = (3' a 3 /(l + al) 3 , (1.5a) 
C = f3lal/(l + alf . (1.5b) 

We see explicitly in Eqs. ( |1.4j ), ( |1.5| ) that all deviations from general relativity tend to zero 
with a a ^ least as fast as Qq. This holds true for weak- field deviations of arbitrary post- 
Newtonian order J7|. Therefore, light-deflection or time-delay experiments || which set 
(through Eq. ( |1.4a|) ) the following limit on the coupling strength of the scalar field, 

al < 1(T 3 , (1.6) 

tightly constrain the theoretically expectable^ level of deviation from general relativity in 
all other experiments probing weak gravitational fields. Note that, in many physically 
motivated models, there are much tighter limits on al coming from equivalence principle tests 
(see e.g. [f| which gets al ^ 10~ 7 in string-derived models). These improved limits crucially 
depend, however, on the detailed structure and magnitude of equivalence-principle- violating 
effects (and disappear in the sub-class of metrically coupled theories). To stay model- 
independent, we shall use the post-Newtonian-derived limit ( |1.6| ) as our standard weak- field 
limit. As we shall see later, the importance of the nonperturbative effects discussed here is 
not uniformly decreased when ao takes smaller values, but can level off or even be amplified. 



In a previous work [10], we have shown that experiments involving the strong gravitational 
fields of neutron stars can exhibit a remarkably different behavior from weak-field solar- 
system experiments. We proved that when a certain mild inequality restricting the curvature 
of the coupling function \nA((p) was satisfied, namely 



We assume here the absence of unnaturally large dimensionless numbers appearing in the suc- 
cessive derivatives of a(ip) : Pq, (3' , . . . 
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A ^ S -4 , (1.7) 



nonperturbative strong-gravitational-field effects developed in neutron stars and induced 
order-of-unity deviations from general relativity, even for arbitrary small values of the linear 
coupling strength Qq. The aim of the present paper is to further study these nonperturbative 
phenomena and to prepare the ground for a systematic application to binary-pulsar exper- 
iments [|TTJ by computing the observational effects depending upon the inertia moments of 
neutron stars. One of the main results of the present study will be to show explicitly that 
binary-pulsar experiments are, in some regions of theory space, much more constraining 
than solar-system experiments. This will be illustrated in an exclusion plot discussed below. 
The organization of this paper is as follows. In Section II, we show how the non- 



perturbative scalar-field effects discovered in JTTJ can be interpreted as a "spontaneous 
scalarization" of neutron stars, analogous to the spontaneous magnetization of ferromag- 
nets. We write in Section III the field equations that must be numerically integrated to 
study these non-perturbative effects in slowly rotating neutron stars. Section IV discusses 
the "gravitational form factors" governing the physics of neutron stars in tensor-scalar grav- 
ity, notably a parameter linked to the variation of a pulsar's inertia moment caused by the 
presence of an orbiting companion. The constraints imposed by three binary-pulsar experi- 
ments on a generic class of tensor-scalar models are then derived in Section V. Finally, the 
conclusions of our study are given in Section VI. 



II. SPONTANEOUS SCALARIZATION 

Before tackling the technical problems posed by the computation of various gravitational 
"form factors" in presence of strong-scalar-field effects, let us clarify, at the conceptual level, 



the physical origin of the nonperturbative effect discovered in [10]. 
Let us consider a very simple coupling function of the form 

A(<p)=Ap(<p) = e X p(±/3<p 2> ) , (2.1) 

corresponding to a coupling strength a((p) = d In A((p)/ dip = flip, where (3 is a given pa- 
rameter. The model ( |2.1| ), where In A(<p) is quadratic in (p, is second in simplicity to the 
Jordan-Fierz-Brans-Dicke model where In A(<p) = a^p is linear in p. [We shall sometimes 
refer to ( f2.1| ) as "the quadratic model".] When (3 satisfies (3 < —4, we are in a regime 
where nonperturbative effects develop for massive enough neutron stars. The results of 
lOH raise a paradox in the limit where the asymptotic value of po tends toward zero, i.e., 



a = (3ip — > 0. Indeed, in the case a = (3p> the right-hand side of Eq. ( |1.2b| ) is proportional 
to tp, and <p(x) = is an exact solution which satisfies the homogeneous boundary condi- 
tions <p — > at spatial infinity. Eq. ( 1.2b|) being elliptic in the stationary case of an isolated 



star, it would seem that the solution, with given boundary conditions, must be unique, and 
therefore that in the homogeneous case <pq = the only solution must be the trivial one 
<p(x) = 0. This conclusion is correct in the case of weakly self-gravitating systems (such 
as ordinary stars, white dwarfs or even low-mass neutron stars). Should not then physical 
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continuity require to take always as "correct" solution of Eq. ( |1.2b| ) the trivial one, even 
when considering strongly self-gravitating systems such as neutron stars ? What can cause a 
discontinuity in the configuration of the scalar field (with homogeneous boundary condition) 
for massive neutron stars ? In the simple case of the coupling function fl2.ip , we have the 
further paradox that the theory is symmetric under the reflection (p —>■ —p, so that it seems 
at face value that the solution of Eqs. (|1.2|) corresponding to the self-symmetric boundary 
conditions (po = must be self-symmetric and therefore identically zero. 

A solution of these paradoxes, and a clearer understanding of the phenomena studied 
rUf, is obtained by making an analogy with the well-known phenomenon of spontaneous 



in 



magnetization of ferromagnets (below the Curie temperature). In the latter case, a conve- 
nient order parameter is the total magnetization M (which is thermodynamically conjugate 
to the external magnetic field Bo : M = — dE/dUo). In our "scalarization" case, we can 
take as order parameter the total scalar charge uja developed by the neutron star (labeled 
A) in presence of an external scalar field ipo; it is defined as the coefficient of G*/r in the 
far scalar field around A : ip(r) = p + G^uja/t + 0{l/r 2 ) as r — > oo. As shown in ||, u>a 
is energetically conjugate to the external scalar field tpo, 

u A = -dm A /d(p , (2.2) 

where rriA denotes the total mass-energy of the star (in Einstein units). It is also the 
quantity which appears directly in the Keplerian-order interaction energy between two stars: 
Mnt = —G*mAmB/rAB — G >a> a<^ 'b / ^ 'ab , where the first term comes from the exchange of 
a graviton and the second from the exchange of a scalaron. In the presence of a non-zero 
external ipo, weakly self-gravitating objects develop a scalar charge which is proportional to 
Po in the limit <p — > ("scalar susceptibility"; the analog to the magnetic susceptibility 
M = xB for weak external magnetic fields in absence of spontaneous magnetization). 

Following Landau, we can understand what happens for strongly-self-gravitating objects 
by writing the total energy to be minimized as a function of both the external field and 
the order parameter, uia^a, y?o) = M^a) — wa<Po, and by assuming that the (Legendre 
transform) energy function ^{oJa) develops, when some control parameter varies, a mini- 
mum at a non-zero value of uja- In our case, if we fix the shape of the coupling function 
A{ip) (for instance Eq. fl2.1|) with /3 sufficiently negative), the control parameter is the total 
baryon mass Tua of the star. A simple model exhibiting the appearance of a "spontaneous 
scalarization" of a star in absence of external field ipo is simply the usual Landau ansatz 
near the critical transition point: u,{uja) = \a(m CT — 1tia)uj 2 a + \boj\- m absence of external 
field, po = 0, the energy tua is minimum at the unique (trivial) solution uja = when 
rriA < while when TtTa > "T CT , there appear two energetically favored nontrivial solu- 
tions uja = ±[6~ 1 a(m j 4 — m cr )] 1//2 . At the critical transition Wa = m CI , the slope dujA/drfiA 
is infinite. As in the ferromagnetic case, the presence of an external field (po ^ smoothes 
the transition. For instance, the "scalar susceptibility" xa = duA/dpo which blows up 
near the critical point as \rf%A — ^cr| _1 when ipo = becomes a rapidly varying but smooth 
function of Wla when ipo ^ 0. The results of [ffO] clearly exhibit the sharpening of the tran- 
sition as po — > 0. This is illustrated in Fig. [I], which displays two curves corresponding to 
Po = 2.4 x 10~ 3 and <^ = for the same theory (/? = —6 in Eq. ( |2.1| )) and the same equation 
of state (EOS II of Ref. |12| ). Note that, when <p ^ 0, it is the sign of the external <p which 
determines the direction of the symmetry breaking. 
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It is convenient, notably for the applications to binary-pulsar experiments, to replace the 
quantity uj a by the related quantity 

u A d\nm A 

a A = = — , (2.3) 

m A dtfo 

which measures the effective strength of the coupling between (p and the star. It is the 
strong-field counterpart of the weak-field coupling strength a = a(<po) and reduces to it 
in the case of negligible self-gravity. Correlatively, it is convenient to replace the scalar 
susceptibility \a = du A /dip by the quantity 

Pa = %±, (2.4) 

which is the strong-field analog of the quantity j3 = da((p )/dip entering the Eddington 
parameter /?Edd — 1, Eq. ( |1.4b| ). The quantity [3 A directly enters many observable orbital 



effects in binary-pulsar systems ||. 

Summarizing, we conclude that the nonperturbative phenomenon discussed in [ITJ can 
be simply interpreted as a "spontaneous scalarization" phenomenon, i.e., a scalar analog of 
ferromagnetism. The condition for this phenomenon to occur in actual neutron stars depends 
on the equation of state of neutron matter. For a polytropic model representing a realistic 
equation of state (with maximum baryonic mass of 2.23 m & in general relativity), we found 
that the critical baryonic mass^ for spontaneous scalarization is smaller than about 1.5 m© 
(which corresponds to a general relativistic mass ~ 1.4 m Q ) when /3 Q = d 2 In A((p ) / dip 2 , < 
—5. For such values of /?o, actual neutron stars observed in binary pulsars would develop 
strong scalar charges even in absence of external scalar solicitation (i.e., even if a = a(y?o) = 
0). For values —5 < /3q < —4, one can still obtain important deviations from general 
relativity if the cosmological value of «o saturates the present weak- field limit ( |1.6|) . In all 
cases, the presence of a non-zero external «o smoothes the phase transition and leads to 
continuously (but fast) varying values of the effective coupling parameters a A and (3 A as 
functions of the mass. Fig. || displays the dependence m CY ({3) for the quadratic model fl2.1| ). 
Some representative numerical values are quoted in Table |. For (3q above some critical value 
/3 cr w —4.34, the maximum mass is reached before the zero-mode can develop. It is plausible 
(but difficult to confirm numerically) that as (3 — > (3 cr , the critical baryonic mass tends to 
the general relativistic maximum baryonic mass (~ 2.23 m & in our polytropic model). 

The behavior discussed above concerns the scalar models invariant under the reflection 
symmetry ip — > —<p, such as A((p) = exp(|/3v2 2 ) or A(ip) = cos(v / — Pf)- A dissymmetric 
coupling function, such as A(ip) = exp(|/5y9 2 + |/?V 3 ) would lead to hysteresis phenomena 
(first-order rather than second-order phase transition): for some values of the control param- 
eter m A , there will be two locally stable energy minima available. The scalar configuration 



2 Note that one can determine the critical baryonic mass as a function of 0, in the quadratic model 



(2.1), by solving a linear problem. Indeed, the onset of the transition happens when Eq. (1.2b) 
with a((p) = f3(f (and and T* replaced by a background general relativistic solution) first admits 
a "zero-mode", i.e., a nontrivial homogeneous solution with vanishing boundary conditions |1C]. 
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chosen by the star would depend on the route taken to evolve into its present mass state. 
Let us also mention that we would get an even richer (Goldstone-like) phenomenology if we 
were to consider models involving several scalar fields, with e.g. spontaneous breaking of a 
continuous symmetry in the scalar-field space. Finally, let us make it clear that a negative 
value for /3q = d 2 In A/ dip 2 , does not mean at all that we are introducing some pathology in 
our scalar-field model. The theories we consider are well-behaved field models having only 
positive-energy excitations. A negative value of /3q means only that scalar field nonlineari- 
ties can reinforce the naturally attractive character of scalar interactions, so that it becomes 
energetically favorable to generate more scalar-field energy^. 

III. SLOWLY ROTATING NEUTRON STARS IN TENSOR SCALAR GRAVITY 

One of the main objects of the present paper is to show how to compute the moments 
of inertia of slowly rotating neutron stars in tensor-scalar gravity, especially in presence of 
the nonperturbative strong-scalar- field effects recalled above. We shall work in the Einstein 
conformal frame, within which the basic global mechanical quantities, such as total mass 
and total angular momentum, are conserved (in absence of radiation or particle exchange) 
and can, as usual, be read off the asymptotic expansion of the metric. The total mass 
(in Einstein units) can be read off the 1/r behavior of g^ or g*,, while the ^-component of 
the total angular momentum J a (in Einstein units) can be read off the 1/r 2 behavior of the 
mixed component g^. We consider only stationary axisymmetric field configurations. It has 
been shown by Hartle |13| (see also [[14]]) that the metric corresponding to a slowly rotating 
star could be written, when keeping only first order terms in the angular velocity Q = u'^/u t , 
as 

dsl = g^dx^dx" = -e v{p) c 2 dt 2 + e^ p) dp 2 

+p 2 d6 2 + p 2 sin 2 e{d<j) + [w(p, 9) - tt]dtf . (3.1) 

Thanks to the neglect of fractional corrections of order Q 2 , the diagonal metric coefficients 
v{p) and /i(p) can be taken to be the solutions corresponding to a spherically symmetric non- 
rotating star. The only new field variable which appears in the slowly rotating case is the 
function u(p, 0) entering the mixed component = p 2 sin 2 9[u(p, 6)— Q]. The subtraction of 



3 The appearance of a negative critical value of (3q can be easily understood in the lowest approx- 
imation, where the scalar energy functional to be minimized reads (when setting G = c = 1): 
E[<p] = I d 3 x g^((9j(/?) 2 + p(l + \f3tp 2 ) ■ Indeed, let us consider for instance the simplest trial 
continuous field configurations, (p(r) = const. = uja/R inside a star of mass m = J d?x p, and 
tp(r) = uja/t outside the star (r > R). This yields = m + ^Cur^, where C = R^ 1 (1 + (3m / R) 

becomes negative for a sufficiently negative (3. The missing stabilizing contribution +^buj\ would 
come from taking into account higher-order nonlinearities. 
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the star's angular velocity Q is chosen for later convenience 4 . The total angular momentum 
J a is read off the 1/p 3 behavior of u : 

« = n-^+0(l). (3.2) 



Then the inertia moment (in Einstein units) is defined, in the slow rotation limit, as the 
ratio 

I A = ^ + 0(n 2 ). (3.3) 

We need now to write down explicitly the field equations ( |1.2|) . As the scalar field <p does 
not couple linearly to the rotation, the field equation for <p is, modulo terms of order Q 2 , 
the same as for a spherically symmetric, non- rotating star (therefore <p will, modulo 0(Q 2 ), 
be spherically symmetric). The field equation for the new variable u comes from 

R% = 2d m ?d a p + . (3.4) 

Simply from axisymmetry (d^ = 0) we see that the scalar contribution to the right-hand 
side of ( |3.4j ) vanishes exactly. We are then left with the usual Einstein field equations with a 
localized material source. Taking as usual a perfect fluid description of nuclear matter (with 
energy density e* and pressure in Einstein units) we can directly use the results of Refs. 
T3|]14|| . [One must, however, be careful not to use equations where the "diagonal" Einstein 



field equations have been replaced.] We find the following homogeneous equation for lu : 



7 dp 



i e -(u +l >)/2 d 1 + ^-^d e (sin'Odeuj) 



IC7ra (e,+p,)e^)/ 2 a;. (3.5) 



c 4 



As in Refs. [13|,14|, a decomposition of u(p,6) in associated Legendre polynomials 



dPi(cos 6) I d cos 9 shows that there is only a P contribution (£ = 1), so that, in fact, lu 
depends only on p and not on 9. Adding the scalar-modified diagonal Einstein equations 
(written in [[LUfl), we finally get the following complete set of radial equations for our field 
variables (a prime denoting d/dp): 

M> = ^p?A'(<p)e + ip(p - 2M)i, 2 , (3.6a) 

c 4 p-2M Hy p(p-2M) ' v ; 

p' = ip , (3.6c) 



4 With this definition of variables, the stress-energy tensor of the fluid gives simply T*^ = (e* + 
v*)e~ u p 2 uj sin 2 9 thanks to a combination between gi± and = v$ /u 1 . 
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p 

M 

oo 
w' 



47rG, pA A (if) . . s . ^ 2(p-M) , 

1 ' [a(tp)(e-3p)+fxP(e-p)}- ^ n J> 



c 4 p-2M 



47rq p 2 A%)p 
c 4 p - 2M 



+ oW> + 



M 



p(p-2M) 



p(p - 2M) 



AnrribnA 3 (ip) 



1 - 2M/p 



w , 



/r> A 4 M(e + p)f ro + ^U^ 2 p--W. 



p-2M 



P 



P, 



(3.6d) 
(3.6e) 

(3.6f) 

(3.6g) 
(3.6h) 



The notation used in Eqs. ( |3.6|) is the following: M(p) is defined by writing the radial metric 
coefficient g pp as = (f — 2M(p)/p)~ 1 . As usual the value of M(p) at infinity is the total 
(ADM) mass. The fluid variables have been expressed in physical units using T" = A 4 (<p)T^. 
[It is in these units that one can write a usual equation of state e = e(n), p = pin), 
where n denotes the physical number density of baryons.] if) and w are just intermediate 
notations for the radial derivatives of <p and u, respectively. Finally, we have added an 
equation for the radial distribution of the baryonic mass Wa = M(R) = rht, f A n,yfg u°d 3 x = 
rfib J ( ^4:7rnA 3 ((p)p 2 (l—2M/p)~ 1 / 2 dp J where R denotes the (Schwarzschild-coordinates) radius 
of the star (i.e., the value of p where p and e vanish). 

Note that several of the right-hand sides of Eqs. ( |3.6| ) contain terms proportional to 
ip 2 = (<p') 2 (i.e., proportional to the scalar-field energy density). These terms do not vanish 
outside the star. However, one can avoid numerically integrating Eqs. ( |3.6|) up to p = oo 
by matching the result of integrating (|3.6|) up to the radius R of the star to the known 
general form of the exact static, spherically symmetric exterior solution. This is, however, a 
bit subtle because the general exterior solution can only be written in closed form in some 
special coordinates introduced by Just [15,1G,d|, or (through a simple transformation) in 
isotropic coordinates, but not in the Schwarzschild coordinates we are using. Still, it was 
shown in [ID] how to extract, via a matching across the star's surface, the global quantities 
rriA and a a from the results of integrating Eqs. ( |3.6| ) up to p = R. We need to do here more 
work to extract J a (and I a) from the results for the variables to and w = duj/dp. 

Outside the star, Eq. (|3|) (with d B = 0) shows directly that p*e~ {u+ ^ /2 d p uj is a constant. 
From Eq. ( |3.2|) , this constant is simply related to the total angular momentum, so that 



du G* T e^' 2 

T = 6 ~ Ja i — 

dp c z p 4 



(outside the star). 



(3.7) 



Eq. ( p.7|) gives one equation to determine J a- We need another equation to determine f2 
and then Ia = Ja/O. Note that the equation for u (e.g. Eq. (|3.5|) ) is homogeneous in uj. 
Therefore, we can start the radial integration with an arbitrary (non zero) value of oj(p) at 
p = 0, but we need to extract from u(p) the value of the fluid angular velocity Q implied 
by this arbitrary choice. To achieve this, it suffices to integrate explicitly Eq. ( |3.7| ) with 
the boundary condition oo(p) — > Q when p — > oo (as is clear from Eqs. (|3.1|) or ( p. 2D ). This 
integration can be done by rewriting Eq. ( p. 7|) in Just radial coordinate r. Indeed, the 
general exterior static, spherically symmetric solution [|l5|,[il||| reads 
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dsi = -e u c 2 dt 2 



+e~ 



dr 2 + (r 2 - ar)(d6 2 + sin 2 9dcf> 2 ) 



r 

i \ d i f-, a 

<P{r) = v?o + - In 1 

a \ r 



(3.8a) 
(3.8b) 
(3.8c) 



where the integration constants a, b, d are constrained by a 2 — b 2 = Ad 2 , and are expressible 
in terms of the total Einstein mass and the effective coupling constant a^, Eq. fl2.3| ), 
via 



b = 2—m A 



a 

b = 
d 1 

— = —OLA ■ 

b 2 



(3.9a) 
(3.9b) 

(3.9c) 



Comparing Eq. ( |3.8a|) with the Schwarzschild form ( p.l|) yields 



a \ (a-b)/2a 

p = r [ 1 



1 - 



a 



a + b s 
2r 



(3.10a) 
(3.10b) 



Inserting these results into Eq. (|3.7|) leads to an elementary integral for ui{r). To write 
explicitly the answer it is convenient to introduce the parameter 



1 / a 
p = — In 1 

a \ r 



In terms of p, the exact exterior solution for u reads 

6G*Ja 



(3.11) 



to = n + 



X 



c 2 b(Ab 2 - a 2 ) \ 

2 

(cosh(ap) — 1) 



1 + e 2bp x 



2b" 
a 



2b 



sinh(ap) 



(3.12) 



Combining the results just derived on the radial dependence of to with the results of |10 
for the matching of the other field variables, we can finally write a set of equations allowing 
one to extract all the needed physical quantities from the surface values obtained from 
integrating Eqs. (|3.6|) from the center p = : 



R 



Ps 



o 2M S 
v' s = Ri, 2 s + 



R(R - 2M S ) 



(3.13a) 
(3.13b) 
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a A = -7— (3.13c) 



g! = (l + ay 1/2 , (3.13d) 
g 2 = (l-2M s /i?) 1/2 , (3.13e) 



^^-^ arCtanh (l + 2^-i) ' (3 ' 13f) 



1 ^ . 
(^0 = ^ - t^a^s ' (3-13g) 



2 

-^-m A = -u' s R 2 Q 2 exp(^-V s ^ , (3.13h) 

m A = M s , (3.131) 

-^J A = ^ro s i? 4 Q 2 exp > (3-13j) 

c 4 3J A f 2i > 3 , 4G*m A 



" = " a " G2 4mi (3 - a\) V ~ 1 + ~R^~ * 



i?c 2 



e^ 2 cosh Qqa)]} 



(3.13k) 

- # • (3.131) 



The notation used in Eqs. ( |3.13| ) is that a suffix s denotes the surface value of any of the 



variables entering the first-order system ( |3.6|) . The only exception (apart from v' s that we 
redefine explicitly as the surface value of the right-hand side of Eq. ( |3.6b| )) is 9 S , which is the 
"correct" value of v at the surface when v is normalized as being zero at infinity. Indeed, as 
the system Q3.6p is integrated from the center (starting with an arbitrary value of v(0)) up 
to the surface, the surface value of v{p) naively obtained from integrating (|3.6| ) is not the 
one to be used in any of the physically normalized results. 

Let us finally mention the set of initial conditions, at the center, used for integrating 
Eqs. ( p.6|) . Actually, because of the singular nature of the point p = 0, one numerically 
imposes initial conditions at a small but nonzero radius p m i n . The values of some of the radial 
derivatives (<p' = ip and uj' = w) are determined so as to be consistent with regular Taylor 
expansions at the origin (for instance, writing <p(p) = <p(x) = <fi(0) + |x 2 A<p(0) + 0(x 4 ) 
determines (p'(p) ~ |pA<p(0) as p — > 0). The complete set of initial conditions reads: 

(3.14a) 
(3.14b) 
(3.14c) 



M(p min ) 


= , 


KPmin) 


= , 


<p(p m in) 


= <Pc 


^(Pmin) 


=(b 




xa 


n(p TO in) 


= n c . 


M(p min ) 


= 0, 



4ttG 



x 



c 4 



(3.14d) 
(3.14e) 
(3.14f) 
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^(Pmin) = 1 , (3.14g) 
^(Pmin) = (^-Pmir^j X ^— -A A (<f c ) 

x [e(n c ) + p(n c )] w(p min ) . (3.14h) 



Note that (as discussed above) the initial conditions (|3.14b| ) and (|3.14gj) are arbitrary, and 



that we transform Eq. ( |3.6ej ) in an evolution equation for the physical number density n 
using the equation of state, i.e., p' = (dp(n)/dn) x n' . The choice of (p c and n c is discussed 
below. 



IV. THE GRAVITATIONAL FORM FACTORS OF ROTATING NEUTRON 

STARS 

A. Scalar-field dependence of the inertia moment 

Extending the analysis of ||10||, we have studied the impact of scalar-induced strong-field 



effects on the gravitational form factors of neutron stars. By "gravitational form factor" 
we mean the set of coupling constants that appear, within tensor-scalar theories, in the de- 
scription of the relativistic motion and timing of binary (and isolated) pulsars. As discussed 
in detail in ||[| the (v /c) 2 -accurate orbital dynamics of binary systems depends, besides the 
Einstein masses of the two objects rriA and m^, on the effective scalar coupling constants 
a a, &b, defined in Eq. (|2.3|), as well as on their scalar- field derivatives Pa, Pb, Eq. ( fO[ ). It 
was also shown in || that the same parameters a a, &b, Pa, Pb, suffice to express all radiation 
reaction effects (up to 0(v 7 /c 7 )) in a tensor-scalar description of compact binary systems. 
On the other hand, the relativistic timing of binary-pulsar systems involves, besides the 
above a's and /3's, a new parameter describing the possible field dependence of the inertia 
moment I a of the pulsar. [In the following, we use the label A to indicate the timed pulsar, 
by opposition to the companion labeled B] Indeed, as pointed out by Eardley [[17] 



sec 



also Jl8|]), the adiabatic invariance (under the slow variation of the local scalar- field environ- 
ment caused by the motion of the companion) of the total angular momentum of the pulsar 
J a = Ia^oaJ^a implies that the angular velocity of the pulsar Qa will fluctuate in response 
to the orbital-induced variations of the external scalar field locally felt by the pulsar. 
As discussed in more detail below, the observable deviations from general relativity implied 
by this effect are given by the parameter = —aBd\ialA/dip , in which I a denotes, as 
above, the inertia moment of the pulsar in (local) Einstein units. 

To compute <91n iA/dipo, we have numerically integrated equations fl3.6p with a suitable 
"shooting" strategy for the choice of initial conditions. Indeed, the quantities that are 
physically fixed are y? (the value of <p far from the star) and tua (the baryonic mass of the 
neutron star). [Note that when a derivative with respect to <po is taken, as in the definitions 
of Pa, Eq. ( ^4|) , or of K^, it must be performed for a fixed value of m^.] Therefore, by trial 



5 We restrict here the more general results of (5| to the simple case where there is only one scalar 
field. 
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and error, one must vary the initial conditions tp c and n c in Eqs. ( |3.14|) until they lead to 



the desired values of ^ and Wla- In the end, one wants to explore the way the observables 
rriA, a a, Pa, I a, 91a Ia/9(Pq, . . . depend upon ip and rf%A- 

The values of ttia, a a, ... as functions of <p and depend upon the equation of state 
used to describe the nuclear matter in the neutron star. We shall discuss in a later publication 
the dependence of our results on the choice of the equation of state. In the present work, 
we shall consider, for simplicity, only a fixed polytropic equation of state: 

_ — . Kn m b ( n\ r 

p = Kn m b ( — J • (4.1b) 



All quantities in Eqs. ( |4.1| ) are in local physical units; m b = 1.66 x 10~ 24 g is a fiducial 



baryon mass and no = 0.1 fm 3 a typical nuclear number density. We shall use the following 
specific values of the polytropic parameters T and K, 

T = 2.34 , K = 0.0195 , (4.2) 

which have been chosen to fit a realistic equation of state which is neither too hard nor 
too soft: the equation of state II of Ref. []I2"| . [The polytropic constant K should not be 



confused with the parameter K% linked to the scalar-field-induced variation of the inertia 
moment.] The precise values (f4.2|) were adjusted to fit the curve giving, in general relativity, 
the fractional binding energy / = (m — m)/m as a function of the baryonic mass. In 
particular they lead to the same maximum baryonic mass, m max = 2.23 m , in general 
relativity. Let us note in passing that to convert from the nuclear fiducial quantities to more 
adequate astrophysical units (m for masses, G*m /c 2 for distances), it is convenient to use 
the numerical value 



A-nG*n Q m b /G*m x 1 



c 2 V c 2 J 296.135 



(4.3) 



For technical convenience, when comparing different theories we keep fixed G* = 6.67 x 
10 _8 cm 3 g _1 s -2 (and m = 1.99 x 10 33 g, measured in g* units). See Ref. for the factors 
(differing from unity by ^ 10~ 3 ) relating g*-frame quantities to directly observable ones. 

We present in Figure ^ some of our numerical results for the dependence upon the 
baryonic mass of a a-, Pa, I a [ m units of m Q (G*m & /c 2 ) 2 ] and din Ia/9(p . All the results of 
these Figures have been computed within the tensor-scalar theory defined by the particular 
coupling function 

M<P) = exp(-3v? 2 ) • (4.4) 

This model belongs to the class of quadratic models ( |2.1| ), and possesses a curvature param- 
eter for the logarithm of the coupling function, (3 = (3q = d 2 In Afd^p\ = —6. In the limit 
where tpo — > 0, this model exhibits a spontaneous scalarization above a critical baryonic mass 
m cr = 1.24 m . As explained in Section II, the presence of a nonzero external scalar back- 
ground ipo 7^ smoothes the scalarization and leads to continuous variations of ctA, Pa, ■ ■ ■ 
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in function of tua- For instance, instead of having a Curie-type blow-up oc \rriA — ?7icr| 1 
for the zero-external-field "susceptibility" Pa = daA/dfo, we get a "resonance" bump in 
Pa when Wla ~ rn cv . There remains however an infinite blow-up in Pa when In a reaches 
the maximum baryonic mass. It is easy to see analytically that this blow-up must be there. 
(The same remark applies to dlnlA/dipo.) For defmiteness, we have drawn Fig. for the 
value 

y? = <^ ax = 2-4 x 10" 3 , (4.5) 



which is the maximum value of ipo allowed by present weak- field tests within the model (JOj). 
This maximum value is obtained from considering not only the limit ajj < 10~ 3 , Eq. ( |1.6[ ), 
coming from time-delay and light-deflection experiments ||, but also the limit 

\Po\al < 1.2 x 1(T 3 (4.6) 



coming from the lunar-laser-ranging constraint |/3| < 6 x 10~ 4 ]19|] on the Eddington parame- 
ter P = /?Edd — 1 ~ |/5o«o ( see Fq- ( P-.4b|) ). When \p \ > 1.2, the limit ( (L6|) is more stringent 
than ( |1.6| ) and defines the maximal allowed value for |a | and thereby for \{p \ « \a /p \ (see 
the exclusion plot in section V.D. below). 

Besides the variation of the shapes of the curves in Fig. ^| when </?o is allowed to vary 
(which is always a sharpening of the bumps and a stabilization of the other features^), we 
have also numerically explored the effect of varying the curvature parameter P in Eq. ( p.l|) . 
The two main effects of varying P are (i) to enlarge the values of the form factors \(Xa\, 
\Pa\, \dha. Ia/9(Po\ as —/3 increases, and (ii) to displace the location of the critical point 
m cr . For instance, we find (within the models Q2.1Q ) m CI (P = —5) = 1.56 m , m cr (/3 = 
—4.5) = 1.84 m Q . These values are below the (expected) maximum mass of a neutron star. 
However, observed neutron stars have baryonic masses around 1.5 m Q (corresponding to 
general relativistic Einstein masses around 1.4 m ), therefore we expect that strong-scalar- 
field effects can have significant observational consequences only when P < —5. 



B. Scalar-field effects in the timing parameter 7 



Up to now, the non-Einsteinian effects linked to the field dependence of the inertia 
moment have been treated by an approximation ||i7| , |i~8| ,|5|| which is insufficient for tackling 
the nonperturbative phenomena discussed here. One of the main aims of the present paper 
is to remedy this situation. Let us first clarify the observable effect of the variation of the 
pulsar inertia moment with the local scalar backgroundF] if a 



6 See, for instance, Fig. [| above which shows that the wide plateau in cxa, beyond m cr , varies very 
little when tends to zero. 

7 This denotes the nearly uniform value of ip on a sphere centered on A having a radius much 
larger than the radius of the neutron star A but much smaller than the distance to the companion. 
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The central tool of binary-pulsar experiments is the "timing formula" (see e.g. [f20|,pT|j ) , 
i.e., the mathematical function relating the "intrinsic time" of the pulsar clock T to the 
arrival time on Earth of radio pulses. The successive ticks of the pulsar time T are defined 
to correspond to successive 2tt rotations of the pulsar around itself: PSR = 2nT/P p , where 
P p is the intrinsic period of the pulsar (for simplicity we neglect here the slow-down of the 
rotation of the pulsar as well as aberration effects). In other words, adding the label A 
and passing to a differential formulation, dT A = Cd(j) A for a certain constant C. In (local) 
Einstein units, the pulsar angular momentum reads J a = Ia^a = IaA^a/ drX-, where dr A = 
\ds* A \/c = (-g*^dz^dz u A ) 1/2 /c is the Einstein proper time in a local inertial frame around 
A. The angular momentum J a is an action variable (J — — -^§Pidq 1 ) and therefore 
an adiabatic invariant under slow changes of parameters. It remains therefore constant as 
the pulsar moves on its orbit and feels a slowly changing if a from its companion. This 
yields gO~a = C'dr A / Ia for some new constant C. The latter equation can be approximately 
rewritten in terms of some coordinate time t used to describe the binary motion: 



dT A « C'V^V 1 - v 2 a/c 2 dt/I A (Mt)) , (4.7) 

where (to sufficient accuracy) v A is the Euclidean square of the coordinate velocity of the 
pulsar = dzA/dt. Using (see ||) 



^ -^ +o (i), „ 8b) 

and the standard relations given by Newtonian orbital dynamics (with effective Newtonian 
constant Gab = G*(l + «a«b))j we find a usual "Einstein" contribution, A E = 7sinii, 



to the timing formula In A#, u denotes the function of T4 defined by solving 



u 



esinw = 2vr[(T A - T )/P b - \P b ((T A - T )/P b ) 2 ], and| 
7 = 7 th (m A ,m B ) 

\ 2/3 

e X B ( G A B(m A + m B )n\ 



n 1 + a ao. b \ c 3 J 

x [X B (l + a A a B ) + l+K A r \ . (4.9) 

The timing parameter 7 should not be confused with the Eddington parameter 7Edd- Here e 
is the orbital eccentricity, n = 2n / P b the orbital circular frequency, X B = ms/imA + m B ), 
and the new contribution K A coming from the variation of Ia under the influence of the 
companion B is defined by 

Kl = -«b^ • (4.10) 



8 The notation 7 th (myi, m B ) in Eq. ( [4.9[ ) refers to the theoretical prediction, within tensor-scalar 
models, giving the phenomenological timing parameter 7 as a function of the masses. See below. 
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Note the dissymmetric roles of the labels A and B. It is important, for applications, to rec- 
ognize that the dependence of the correction K% upon the two masses m^m^ is factorized 
(in the single scalar case that we consider here). Accordingly, it might be convenient to 
define the quantity 

k A {m A ) = -dlnI A /d<p , (4.11) 

so that K A ^{r(iA- l rn B ) = /^(m^a^m^). 

The reasoning above (based on the use of the Einstein conformal frame) could be done 
using the "physical" (or Jordan-Fierz) conformal frame. Indeed, the angular momentum is 
independent of the conformal frame (being an action variable). This means I a^a — 1 a^a so 
that the pulsar intrinsic time (which is a conformal invariant, being proportional to the angle 
(Pa) can be equivalently written as^ dT A = C'cIt a /Ia = Calf a/ I a- The calculation is (as 
always) slightly more complicated in the Jordan-Fierz frame and leads to a correction K A 
instead of the K A in Eq. Q4.9|) , given by the sum of two terms: K A = olqolb — a^lni^/c^o, 
which is (as it should) found to be identically equal to K A , Eq. ( |4 . 1 0| ) , when using the link 
T A = A((p A )I A . 

In previous works [ ]TT| , |i"8| ,|5[| one had assumed, as an approximation, that I a was simply 
a function of the local, externally imposed, value of the effective gravitational "constant" 
G((p A ). Such an assumption is meaningful only in a "quasi- weak- field" approximation where 
one formally considers the compactness of a neutron star as a small expansion parameter 
(see Section 8.1 of ||). This approximation breaks down precisely when the strong-scalar- 
field effects studied here develop (i.e., when > 4). Previous treatments introduced 
the parameter ka = — din Ia/ din G a- When it is meaningfully defined, the parameter ka 
is linked to the parameter k A = —d In I A /dipo introduced above by: Ica ~ «o + 2«o[l + 
/?o/(l + o^ka- This formula shows that, under the assumptions of previous approximate 
treatments, the correction K A was proportional to the product ctoctB between the weak- 
field scalar coupling a and the possibly strong- field amplified effective coupling as- As a 
is observationally strongly constrained, this led always to small values of K A with nearly 
negligible observational effects. By contrast, the exact result ( |4.10| ) is fully sensitive to 



strong-scalar-field effects taking place both in the pulsar and its companion. To illustrate 
the order of magnitude of possible deviations from the general relativistic prediction^ for 
the timing parameter 7 in systems made of two neutron stars, we plot in Fig. ^ the value 
of K A (corresponding to the cases where = m^) versus rriA within the model Ag(yj), 
Eq. fl4.4|) , and using the same assumptions as in Fig. |3] (notably a maximally allowed value 
of (fo, Eq. (fO)|)). We see on Fig. [| that when uia ~ 1 m & , we get very drastic modifications 
of the general relativistic prediction for 7 (except in a small neighborhood of rriA ~ 1.3 m 
where K A vanishes). In particular, when 1.1 < m^/mQ < 1.2, K A takes largish negative 



9 Note in passing that the pulsar clock ticks neither the Einstein time nor the Jordan-Fierz one. 
Indeed, both I a and I a fluctuate because of their dependence on tpA(t). 

10 The general relativistic prediction 7 GR (m 7 4, ms) is obtained from Eq. Q4.9| ) by setting aAO-B = 0, 
Gab = G, and Kf = 0. 
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values which change the sign of the predicted 7 th ! [The minimum value of K% in Fig. f| is 
reached for = 1.13 m and equals K^ mm = —3.45, yielding 7^(777^4,77?^) = — 1.277 GR .] 
We computed also for smaller values of the external scalar field <po and found (as usual 
by now) that they cause a sharpening of the "resonance" bump in Fig. |j. For instance, we 
found that K^ mm = —6.68 for = <£>™ ax /10. Paradoxically, smaller values of the weak-field 
coupling ao predict larger values of the modification to the timing parameter 7, though 
concentrated over a smaller range of mass values. This effects is even more spectacular 
for K*[ when % > 77T4 « m cr : in that case, the effective coupling 0.3 tends to a non- 
vanishing constant as (po — > 0, while dlnl^/difo blows up, so that \K^\ can take arbitrarily 
large values. For instance, one gets K% min = -8.20 for yj = y^ ax , and K% min = -23.82 
for (fo = <^™ ax /10 [yielding 7min( m ^) m s) ~ — 23 7 GR !]. Finally, varying the value of the 
curvature parameter (3 in the models ( |2.1| ) displaces the center of the bump in K^, towards 
lower (higher) values when —f3 increases (decreases). 



V. APPLICATION TO BINARY-PULSAR EXPERIMENTS 

We have now all the necessary tools in hand for exploring the impact of non-perturbative 
scalar effects on binary-pulsar experiments. In a future publication [jll] J, we shall confront 
in a systematic manner the predictions of tensor-scalar gravitation theories with a more 
complete and updated set of binary-pulsar data. In the present work, we shall illustrate how 
binary-pulsar data give us very significant constraints on the strong-field regime of relativistic 
gravity by comparing published data on PSR 1913+16, PSR 1534+12 and PSRQ 0655+64 
with the predictions of tensor-scalar theories exhibiting the nonperturbative effects discussed 
above. In Ref. pT| , we shall also take into account data on certain nearly circular binary 



systems which test the strong equivalence principle p2Tf2~4]j. We do not consider them here 
because they are less constraining than the systems we study. [Indeed, the "Stark effect" is 
proportional to the product a («A — «b) an d, therefore, is already significantly constrained 
by the solar-system limits on ao-] 

The case of PSR 1913+16 is the richest in that it involves many different types of strong- 
field effects: (i) modifications of the first post-Keplerian orbital motion (observable through 
the periastron advance u), (ii) modification of gravitational radiation damping (observable 
through the orbital period decay P&), and (iii) sensitivity of the pulsar inertia moment to 
an external scalar field (observable through the timing parameter 7). As we shall discuss 
below, the case of PSR 1534+12 very usefully complements that of 1913+16 in trading the 
Pb measurement against a measurement of the shape parameter s of the gravitational delay. 
The scalar-field effects in u, Pf, and s have been already worked out with sufficient accuracy 
in the literature P,^,[L8|], while the scalar- field effects in 7 have been discussed above. 



Note that the pulsar community now uses an updated notation in which these pulsars are 
called PSR B 1913+16, PSR B 1534+12 and PSR B 0655+64, respectively. Here, the label B (for 
Besselian) refers to the equatorial coordinate system based on the 1950 equinox [while the letter J 
(for Julian) refers to the 2000 equinox]. 
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A. The PSR 1913+16 experiment 



We recall that, at present, one can phenomeno logically extract from the raw data of 
the binary pulsar PSR 1913+16 (following the methodology of pi] ) three well-measured^ 
observables: uj obs , 7° bs and P fe obs . Here ci; obs denotes the secular rate of advance of the 
periastron, 7° bs denotes the observed value of the timing parameter discussed above, and 
P b ohs denotes the secular change of the orbital period. The values we shall take for these 
observed parameters are |p6 j: 



u ohs 
7 obs 



>obs 



4.226621(11) ° yr~ 
4.295(2) x 1(T 3 s , 
-2.422(6) x 1(T 12 



(5.1a) 
(5.1b) 
(5.1c) 



where figures in parentheses represent la uncertainties in the last quoted digits, 
also need the Keplerian parameters 



We shall 



P b = 27906.9807804(6) s 
e = 0.6171308(4) . 



(5.2a) 
(5.2b) 



The important point (which is the basis of the parametrized post-Keplerian approach 
PH ) is that the observables ( |5.1| ) have been extracted from the raw pulsar data without as- 
suming any specific gravitation theory (at least within the very wide class of boost-invariant 
theories). One can then use the three pieces of data ( |5.1| ) to constrain theories of gravitation. 
To do this one must compute, within the theory to be tested, what are the predictions it 
makes for tit, 7 and P b as functions of the two (a priori unknown) masses m A , m B . We have 
written in Eq. (|4.9| ) above the theoretical prediction for the timing parameter, 7 th (myt, tub), 
within tensor-scalar gravity models. The theoretical prediction for the periastron advance 
rate has been worked out in Refs. [|18| , pl~| ,[5f and reads, with the notation of the present paper, 



u th (m A ,m B ) 



3n 



' G AB {m A + m B )n 



2/3 



1 - \a A a B X A j3 B a 2 A + X B j3 A a B 



1 + a A a B 



6(1 + a A a B ) 



(5.3) 



The notation in Eq. (|5.3| ) is the same as in Eq. ( fOl) . We recall that n = 2n/P b , G AB = 
G*(l + a A a B ), X A = m A /{m A + m B ) = 1 — X B . Finally, the theoretical prediction for the 
(radiation-reaction driven) orbital period decay has been derived in Ref. || (with the full 
needed accuracy and for generic tensor-scalar theories). It is given as a sum of contributions: 



P b ttl (m A ,m B ) 



^imonopolc _|_ ^jdipole _|_ ^jquadrupole 



^pquadrupole _|_ 'pgal 



GR 



+ 5 th ( P 



(5.4) 



12r 



Two more observables, r obs , s obs , are measured with low precision [25]. 
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The first three contributions correspond to energy lost to scalar waves (of monopolar, dipolar 
and quadrupolar type, respectively). The fourth one corresponds to the energy lost to 
quadrupolar tensor waves (pure spin-2 field gt u ). The explicit expressions of these four terms 
are given in Eqs. (6.52) of ||. The fifth contribution is the value of the galactic contribution 
to the observable Pt> computed in within the assumption (true in general relativity) 
that neutron stars fall like ordinary bodies in the gravitational field of the Galaxy. Finally, 
the sixth and last contribution is the modification of the galactic contribution due to the 
fact that, in tensor-scalar gravity, neutron stars fall differently from weakly self-gravitating 
bodies (Eq. (9.22) of §). 

As usual, given a specific tensor-scalar theory, a value for the externally imposed asymp- 
totic boundary condition ip Q , and a specific nuclear equation of state, the three equations 
iu th (mA, tub) = uj ohs , 7 th (myi,mB) = 7° bs , P t th (mA, tub) = Pb hs , define three curves (in fact 
three strips) in the two dimensional plane of the masses {itia^b)- If the three strips meet 
in a small region, the considered tensor-scalar theory is consistent with the binary-pulsar 
data. If they do not meet, the considered theory is inconsistent with the pulsar observations. 

Before presenting the results of such a confrontation for scalar models exhibiting non- 
perturbative effects, it is instructive, for making a contrast, to discuss two other cases 
(besides the general relativistic one): (i) the case of the well-known Jordan-Fierz-Brans- 
Dicke (JFBD) theory and (ii) the case of scalar-tensor models with positive curvature 

of the coupling function. The JFBD theory contains only one free parameter cto and is 
defined by the coupling function 

AjfbdM = exp(a y) • (5.5) 

The scalar coupling strength in this theory is constant: ot((p) = din A/dip = ao- As a 
consequence (see Section II), it cannot exhibit nonperturbative effects. All deviations from 
general relativity, be they in weak-field or strong- field conditions, are proportional to a^, and 
are uniformly constrained by the solar-system limit ( p..6|) . On the other hand, as discussed 
in |IJ, scalar-tensor models belonging to the quadratic class ( |2.1|) with /3 > exhibit 
nonperturbative effects of a deamplification type: deviations from general relativity are 
exponentially suppressed by strong-field effects, i.e., are proportional to a% exp(— 3/?sa), 
where sa ~ 0.2 is a measure of the strength of the self-gravity of the considered neutron 
star. 

As one of the aims of the present work is to delineate the cases where binary-pulsar 
data give more stringent constraints than solar system data on tensor-scalar theories, we 
shall draw the figures below (except when otherwise indicated) under the assumption that 
the externally imposed tpo always takes the maximum value allowed by solar-system tests. 
As said above, the corresponding value of <p depends on the theory considered and is 
determined from combining the two inequalities ( |1.6| ) and ( |4.6| ). Fig. [| exhibits the curves 
defined by the observables u, 7 and Pb in PSR 1913+16 when interpreted in the framework 
of three different theories: general relativity, Jordan-Fierz-Brans-Dicke and the quadratic 



model (|2.1| ) with (3 = +6. This plot illustrates the fact that binary pulsars involving nearly 
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equal mass0 neutron stars do not constrain more severely than solar-system data theories 
with logarithmic coupling functions In A(ip) which are either linear (or nearly linear) in tp 
p8| or have a positive curvature (i.e., convex functions In A(tp)). Note that the Pf, curve 
of a quadratic model with positive curvature lies between the general relativistic and the 
JFBD ones. The fact that in Fig. [5] the P^~ +6 curve almost overlaps with the Co curve is a 
numerical coincidence caused by our choices for (3 and (fo- 

By contrast, tensor-scalar theories involving sufficiently concave functions lnA(ip), i.e., 
models ( |2.1| ) with (3 sufficiently negative, show a very different behavior when confronted 
to pulsar data. This is illustrated in the remaining figures. The left panel of Fig. |B| shows 
that although the quadratic model A-^((p) (i.e., (3 = —4.5 in Eq. ([2.1|)) can pass all present 
solar-system tests, it fails the (uj-j-Pb) 1913+16 test. For such a model, pulsar observations 
constrain more strongly the theory than weak-field tests. This is further illustrated in the 
right panel of this figure, which shows that one needs a smaller value of tpo, i.e., a smaller 
value of «o = ft^Po than the maximal one allowed by solar data. We did not make a 
exhaustive numerical search but it seems that a value ag SR ~ ^ a osofar * s the correct order 
of magnitude that pulsar data can tolerate. Note that, in terms of the basic Eddington 
parameters 7 = 7Edd — 1 ~ — 2«q, (3 = /^Edd — 1 ~ \(3qoI, this means that binary-pulsar 
data are one thousand times more constraining than solar-system tests, constraining 7 and 
(3 (for the considered model Ajj(ip)) below the 10~ 6 level. 

In the above case ((3 = —4.5), the fact that the maximal weak-field-allowed value of 
ipo was forbidden, while a 30 times smaller one was allowed, was due to the presence of 
significant nonperturbative scalar effects for ip = ^[Sofar which tended to zero when <p — > 0. 
This smooth disappearance of nonperturbative effects when ipo — > was connected to the 
fact that the critical value of the mass m cr (above which spontaneous scalarization occurs in 
zero-external- field conditions) is slightly larger that the actual mass of the stars \m cr (f3 = 
—4.5) « 1.84 m G while m^ stRt « m B estRt « 1.5 m Q ]. 

A quite different situation arises for values of —f3 large enough to make m cr smaller than 
the masses of the stars. This case is illustrated in Fig. [7] where one confronts pulsar data 
with the model A§, Eq. ([14]) . [In this model m cr (P = —6) ~ 1.24 m .] This figure shows 
that the model A$ fails very badly the (u-^-Pb) igis + w test, while it can pass all present 
solar-system tests. Especially remarkable is the wild behavior of the 7 curve, which is due 
to the large values of the K\ deviation discussed in the previous section (see Fig. [| there) . In 
that case, this disagreement between theory and experiment is not alleviated by considering 
smaller values of ipo, as illustrated on the right panel of Fig. [I]. From our numerical results, 
we find that the (^-7-A)i9i3+i6 test can be passed, if at all, only for extremely fine-tuned 
values of the masses in close neighborhoods of the critical values w TFTb ~ m cr . Barring 
any fine-tuned coincidence, we conclude that the tensor-scalar theory defined by Ag(ip) is 
incompatible with pulsar data whatever be the value of tpo, even infinitely smaller than y^ofar- 



13 As emphasized by Eardley [{17|| , the situation is different for unequal mass systems thanks to the 
presence of scalar dipole radiation oc (a^ — a_e) 2 . This shows up in Fig. || as a strong distortion of 
the Pb curve away from the diagonal txia = tub- See below our study of the unequal mass system 
PSR 0655+64. 
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This remarkable conclusion proves explicitly that pulsar data are qualitatively different from 
solar-system data in their probing power of relativistic gravity. The theory defined by A$((p) 
could always be made compatible with solar-system tests of any precision^, while it is 
already falsified by existing pulsar observations. As the critical mass m cr decreases when 
—f3 increases, the confrontation theory/pulsar-experiments can only get worse when —(3 > 6. 
Furthermore, our (partial) numerical exploration of the range —6 < (3 < —4.5 shows that 
values of (3 smaller than —5 are already excluded. [This will be illustrated by an exclusion 
plot discussed below.] We conclude that present PSR 1913+16 data already exclude all 
"quadratic" models ( J2.1[ ) with (3 < — 5. 



B. The PSR 1534+12 experiment 

One conceivable deficiency in the above argument is the possible presence of a cosmo- 
logical variation of (po- Indeed, a nonzero value of (po = d(po/dto entails a secular variation 
of the strength of scalar gravity and produces an additional contribution in P b th . From 
the observed P£ hs , one cannot decorrelate 0o-effects from scalar modifications to radiation 
damping. One way to decorrelate 0o-effects is to consider pulsar experiments which are 
insensitive to (po- Such is the case of the measurements of the three observables cl> obs , 7° bs 
and s obs in PSR 1534+12. Here s obs denotes a phenomenological parameter measuring the 
shape of the gravitational time delay f2H]2T|| . The values we shall take for these threeQ 



observable parameters are [23 



u obs = 1.75573(4) yr" 1 (5.6a) 

7 obs = 2.081(16) x 1(T 3 , (5.6b) 

s obs = 0.981(8) . (5.6c) 

We shall also need the Keplerian observables 

P b = 36351.70267(2) s , (5.7a) 

e = 0.2736771(4) , (5.7b) 

x = 3.729458(2) s . (5.7c) 

The theoretical predictions for uj th and 7 th have been written above. That for s th reads PJ 

!lX -[G AB (m A + m B )n/c z }- l l\ (5.8) 



X 



B 



where we have used the same notation as in Eqs. (|4.9|) and ( |5.3| ) above, and where x = a^sjc 
is the projection of the semi-major axis (ai) of the pulsar orbit on the line of sight (in 
light-seconds). 



14 For the argument's sake we assume here that general relativity is the correct theory of gravity. 

3 We do not consi 
fractional precision. 



15 We do not consider here the other observables r obs and P^ bs which are measured with low 
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We have plotted the three curves defined by this (^-7-3)1534+12 test for various values 
of j3 and <p . For instance, we exhibit the case (3 = —6, y? = v^osofar m Fig. || (together 
with the case of general relativity). From our (partial) numerical study, we conclude that 
the quadratic models A@ fail the (^-7-3)1534+12 test when (3 < —5.5. The corresponding 
exclusion plot is very similar to that defined by PSR 1913+16 (see below). 



C. The PSR 0655+64 experiment 

The binary pulsar PSR 0655+64 is composed of a neutron star of mass ~ 1.4 m , and 
a white dwarf companion of mass ~ 0.8 m . They move around each other on a nearly 
circular orbit in a period of about one day. In tensor-scalar gravity, such a dissymmetrical 
system is a powerful emitter of dipolar scalar waves, especially in presence of nonperturbative 
scalar effects. The theoretical prediction || for the corresponding orbital period decay is 
dominated by the 0(v 3 /c 3 ) dipole contribution in Eq. ( |5.4| ) above: 

P?(m A ,m R )~P^ ^m A m B n 



1/3 {rriA + ms)c 3 



1 _|_ e 2 /2 fv 5 \ 

X (1 _ e 2)5/2 - + {^J ■ (5 ' 9) 

The fact that the observed value of P& in PSR 0655+64 is very small (and, in fact, consistent 
with zero) constrains very much the magnitude of the effective coupling strength a A , and 
therefore the possibility of nonperturbative effects. The experimental data we need for our 
analysis are taken from Ref. p3f : 



P b = 88877.06194(4) s , (5.10a) 

e < 3 x 10~ 5 , (5.10b) 

x = (ai sini)/c = 4.12560(2) s , (5.10c) 

P b = (1 ±4) x 10~ 13 . (5.10d) 

The masses of the pulsar and its companion are not known independently. From the observed 
mass function, the a priori statistics of the inclination angle i, and the observed small 
statistical spread of neutron star masses around 1.35 m®, one can deduce a range of probable 
values for the pair (rriA,iTLB) '■ Essentially, one is limited to a sub- region of the rectangle 
rriA = (1.35±O.O5)m , = (O.8±O.l)m in the mass planeP] |p3[ . In our calculations, we 
will choose the mass pair which gives the most conservative bounds on tensor-scalar gravity, 
namely uia = 1.30 m , m# = 0.7 m . 

Finally, using the fact that the self-gravity of the white dwarf companion is negligible 
compared to that of the pulsar (so that ~ «o), we get from Eqs. ( |5.9|) and ( |5.10|) the 1-a 
level constraint 



16 We use here the fact that the scalar modifications to the link between the observed mass function 
n 2 (aisini) 3 and the Einstein masses rriA, ms due to the factor Gab/G = 1 + a^a^ are small 
because as ~ «o for the white dwarf companion. 
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(a A (rn A ) - a ) 2 < 3 x 10 4 



(5.11) 



D. Exclusion plots within a generic two-dimensional plane of tensor scalar theories 

It is instructive to contrast the pulsar constraints on tensor-scalar gravity with the 
constraints obtained from solar-system experiments. We can use the class of quadratic 
models (|2.1|) as a generic description of the shape of the coupling function around the 
current cosmological value of (p. In other words, we can parametrize an interesting class 
of tensor-scalar models by two parameters^} say, cto = «(y?o) and (3q = da(<po) /d(fo. [In 
quadratic models, a = /3ip , and @ = [3 is field independent.] We can then interpret 
all experimental data (solar-system and pulsar ones) as constraints in the two-dimensional 
theory plane (ctojA))- For instance, neglecting the correlations in the measurements of 
the two Eddington parameters 7Edd and /^Edd, solar-system data rule out the regions of 
the (a 0) A)) plane where the inequalities |7Edd — 1| < 2 x 1CT 3 {i.e., < 10~ 3 ) ||, and 
l/^Edd — 1| < 6 x 10 -4 (i.e., \Po\aQ < 1.2 x 1CT 3 ) |19| are not satisfied. The corresponding 
exclusion plot is represented in Fig. []. On the same plot, we can represent the constraints 
brought by pulsar data on tensor-scalar models. The PSRs 1913+16 and 1534+12 data 
give constraints which are numerically similar. Taken together, they exclude the region to 
the left of the wavy line indicated on Fig. || This line is approximate and was obtained by 
interpolating between a few values of cto and (3q. We leave to future work a detailed study 
of the precise region of the (ao, A))-plane excluded by these pulsar data. 

As for the PSR 0655+64 data, they define through Eq. (|5.11| ) (with the conservative value 
itla = 1.30 m & ) another limit on tensor-scalar models. We have numerically computed the 
region of the (a , /?o)-plane defined by the inequality ( |5.11| ). The corresponding allowed 
region is contained within the two dashed lines labeled 0655+64 on Fig. £5[ The region 
simultaneously allowed by all the tests is shaded. 

To prevent any confusion, let us note that the limit on the 2PN parameter ( = (3qO^/ (1 + 
«o) 3 obtained in a recent simplified analytical study of combined pulsar data is valid 
only for \f3\ < 1. Indeed, the approximate treatment of assumed the absence of any 
nonperturbative effect, i.e., an absolute value of (3q appreciably smaller than 4. The PSR 
0655+64 constraint studied here should merge with the limit /?o a o < 4 x 10~ 3 derived in 
when \(3q\ < 1. [We use here the inequality (5.24c) of Ref. [0] together with the theoretical 
constraint that /?o a o be positive.] This merging occurs anyway in a region which is already 
excluded by solar-system data. 



17 This two-parameter class of models is representative of the large class of coupling functions A((p) 
which admit a local minimum and contain no large dimensionless parameters [i.e., we assume that 
higher derivatives f3' = d(3((po)/d(po, (3q = d{3' (ipo)/dipo are of order unity]. 
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VI. CONCLUSIONS 



Before summarizing the main results of the present work, we wish to indicate briefly 
the frameworks within which our findings might be physically relevant. First, let us note 
that (barring any unnatural fine tuning) they are not relevant in the case of a strictly 
massless scalar field ip, or at least when its mass m v <C TiHq/c 2 ~ Ttioo x 2.13 x 10 -33 eV, 
where H = h 100 x 100 km.s _1 /Mpc denotes Hubble's constant. Indeed, in such a case it 



was found, both in traditional equivalence-principle-respecting tensor-scalar theories [29 
and in string-inspired models with a massless dilaton or modulus @, that the cosmological 
evolution naturally drives the vacuum expectation value of ip toward minima of \nA(<p). 
As the latter vacuum expectation value coincides (modulo small fractional corrections due 
to spatial fluctuations of the gravitational potential |7j]) with our externally-imposed <po, a 
natural prediction of these massless models is that oq = d In A((po) / dpo is small and that 
Po — d 2 In A(ipo)/d(p\ is positive. The former result concerning ciq is in agreement with 
observational data, but the latter one concerning j3 gives the wrong sign for spontaneous 
scalarization effects^. 

On the other hand, our results are relevant to the wide class of models comprising scalar 
fields of ranged 10 6 km < A = h/m^c ^ cHq 1 . Assuming that the endemic "Polonyi 
problem" (too much energy stored in the cosmological oscillations of <po(i)) of such 
models is somehow solved |JT| or fine-tuned to provide fl = ptot/Pcr = 1 |32| , the condition 



for these models to be naturally compatible with solar-system constraints is simply that the 
location (p of the minimum of the <p potential, V(<p) = (c 4 /87rG*)m^(<yj — y?o) 2 , coincide 
(or nearly coincide) with an extremum of the coupling function A(<p). Such a coincidence 
would, for instance, naturally follow from a discrete symmetry about (po, say a reflection 
symmetry {ip — (p ) — > —(<p — (p ) |J. Under such a condition, the present value of the 
weak-field scalar coupling a = d\nA((p )/d<po would be naturally extremely small, and the 
sign of the curvature (3q = d 2 In A((p ) / 'd<p\ could be expected to be negative with a priori 
50 % probability. 

As a simple example of such models, we can consider a finite-range scalar field coupled 
only to the gravitational sector through a multiplicative coupling to the scalar curvature, 
say 

167rG* J c v 

+F($)R) + S m [il> m -,g ia ,] , (6.1) 



18 Note, however, that a coupling function of the type In A(<p) = +e 2 p 2 — X 2 p 4 , with a sufficiently 
small e and a sufficiently large A, would reconcile a (very localized) minimum at ip = with a 
mainly concave coupling function leading to nonperturbative effects. We do not wish to consider 
here such fine-tuned cases. 

19 The lower bound comes from the requirement that p> be effectively massless on the scale of a 
typical binary pulsar. If it is violated one must correct our formulas by Yukawa exponential factors. 
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where we have introduced the possibility of an arbitrary field-dependent normalization of 
the kinetic term of $. One transforms (|6.1| ) into our canonical form (|1 . 1| ) [completed by a 
(^-potential term V(ip) = (c 4 /167rG*)F _2 ($)f/($)] by defining 



[IV 1 



3 F' 2 ($) 1 Z($) 

4 F 2 ($) + 2 

This corresponds to a coupling function 

A(tp) = F- 1/2 ($) . 



(6.2a) 
(6.2b) 

(6.3) 



The simplest example of such models is a massive scalar having a nonminimal dimensionless 
coupling to curvature: 



16ttG* ' 

+^$ 2 + i?) + 5 m [^ 



(6.4) 



This corresponds to Z(&) = 1, U(§) = m|$ 2 , and -F($) = 1 + £$ 2 . In that case, one can 
integrate Eq. ( |6.2b| ) explicitly. Assuming for instance that £(1 + 6£) > and introducing 

the notation x = Vs(l + 6£)> one gets 



2^2 (<p-<p ) = - In 



A 



+V/6 In 



1 - 2^^$ 



1 + 2 X $ (yi + x 2 $ 2 + x$ 

yi + X 2 ^ 2 - 
l + £$ 2 



(6.5) 



Note that - <p ) = + 0($ 3 ) when $ -»• 0. From Eqs. dJHEED, one S ets easily 

the following parametric representations 



A(<p)= (l + ^ 2 
din A((p) 



-1/2 



a(ip) 



dip 

-y/2£$ [l + e(l + 60$ 2 
d 2 In A((p) 
dip 2 



-1/2 



-2£ ! + 



x 



l + e(l + 60$ 2 

In particular the value of the curvature parameter around ipo reads 

d 2 In AO) 



A) = /%o) 



-2f ■ 



(6.6a) 



(6.6b) 



(6.6c) 



(6.7) 
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Therefore, positive values of £ [which, in the formulation ( |6.4|) , seem preferred because unable 
to cause a change of sign of the coefficient (1 + £$ 2 ) of the kinetic term for g^ u ] correspond 
to the "interesting" case where spontaneous scalarization effects can occurQ. 
Let us now summarize the main results of the present work: 

• We have clarified the physical meaning of the nonperturbative strong-gravitational- 
field effects discovered in |10| by interpreting them, by analogy with ferromagnetism, as a 



phenomenon of "spontaneous scalarization" . Negative values of the curvature parameter (3q 
(like negative values of the coupling between magnetic moments = g^fij) favor the 
spontaneous creation of a scalar field when considered in the context of gravitationally com- 
pact objects (neutron stars). The critical baryonic mass for the spontaneous scalarization 
transition is ^5 1.5 m when /3 < —5. See Table H and Fig. || for precise values. The pres- 
ence of some externally-imposed scalar field background ipo with «o = d\n A(ip Q )/dipo ^ 
smoothes the scalarization transition (analogously to the presence of an external magnetic 
field for the ferromagnetic transition). 

• The development (through the previous mechanism) of strong scalar fields in neutron 
stars leads to very significant deviations from general relativity. These deviations are mea- 
sured by some gravitational form factors oca, Pa, K a which enter the effects observable in 
binary-pulsar experiments. In our previous work Jill] we focussed on the effective scalar 
coupling constant a a- Here, we gave results for [3a = daA/dip (see also @), and for 



K^y = —asd In Ia/O^q which enters the parametrized post-Keplerian timing parameter 7. 
To compute K^, we generalized to a tensor-scalar context the work of Hartle on the 



general relativistic inertia moments I a of slowly rotating neutron stars. We found that K% 
could cause very drastic deviations from general relativity in tensor-scalar theories contain- 
ing no large dimensionless parameters. This is achieved without fine-tuning and in theories 
having only positive-energy excitations. 

• We presented a preliminary investigation of the confrontation between scalar models 
exhibiting nonperturbative effects and actual binary-pulsar experiments. We contrasted the 
probing power of pulsar experiments to that of solar-system ones by working in the two- 
dimensional (oto = d\iaA({p )/d(p , j3 = d 2 In A(<p ) / <9</9o)-plane describing a generic class 
of tensor-scalar models. Using published data on PSRs 1913+16, 1534+12 and 0655+64 
(and a specific nuclear equation of state), we found that binary-pulsar experiments exclude 
a large domain of theories compatible with solar-system experiments (see the exclusion plot, 
Fig. §). In particular, they constrain /3 (independently of a>o) to 

fa > -5 . (6.8) 

Interestingly, this bound can be expressed in terms of the well-known weak-field Eddington 
parameters 



2o In 

our conventions, the so-called "conformal coupling" corresponds to £ = — g, and to 
a coupling function A(ip) = cosh(tp / y/3) . Note, however, that only the action 5 con f = 
/ d 4 xg 1//2 [— g^d^dpQ — ^R& 2 ] (without bare Einstein term) exhibits a conformal invariance 
p3| . The action S^onf involves a spin-2 but no spin-0 excitation. 
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< 1.3 . (6.9) 

7Edd - 1 

The singular (0/0) nature of the ratio on the left-hand side vividly expresses the fact that 
such a conclusion could never be obtained in weak-field experiments (at least until they find 
a significant deviation from general relativity). It must be kept in mind that the inequality 
( |BT9D is one sided only, and that /?Edd — 1 and 7Edd — 1 must be taken with their signs. 

Let us note that, on the whole, the fact that pulsar data tend to exclude sufficiently 
negative values of /?o is nicely compatible with the expectation from cosmological attractor 
scenarios p9|,|6[| that <p be dynamically driven toward a minimum of In A((p). Though our 
results have been derived by assuming a particularly simple coupling function between the 
scalar field and matter ("quadratic model": In A(<p) = |/?o<^ 2 ), we think they hold in general 
classes of coupling functions containing no small or large dimensionless parameters. In a 
future publication [|ll|], we shall present a more systematic confrontation between tensor- 
scalar theories and binary-pulsar experiments. 

Before ending this paper, we would like to stress some of the limitations of our treatment 
that we intend to overcome in future work: (i) We considered here only one specific equation 
of state, modeled as a simple polytrope; a more complete study should consider a selection 
of realistic equations of state, (ii) In the present paper, we did not consider the effect of a 
cosmological variation of <pq. We are aware that a nonzero 0o of order the Hubble parameter 
could significantly modify the interpretation of some of the pulsar tests discussed above. 
However, as one disposes of several independent pulsar tests, some of which do not involve 
0o-sensitive observables (such as the (u-j-s) 1534+12 test considered above and several "Stark" 



tests [2"2"|-|23|), we think that a combined discussion of all pulsar tests will lead to limits on 
«o and /3q similar to the ones we have obtained here, and, in addition, to significant limits 
on ip . 



Finally, let us note that the strong scalar field effects discussed in |H| and the present 



paper could have a very significant impact on several aspects of the theory of gravitational 
radiation from compact objects (besides the ones taken into account in Section 6.2 of || 
and Eq. (|5.4[) above). Of particular interest would be: (i) the opening of a new, significant 
energy-loss channel in (spherically symmetric!) gravitational collapse and neutron-star bi- 
nary coalescences; (ii) important modifications in the conditions for the onset of radiative 
instabilities in fast-rotating neutron stars (Chandrasekhar-Friedman-Schutz instability) 
Both issues are particularly worth of further study. 
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FIGURES 




FIG. 1. Effective scalar coupling strength — a a = ^a/ m A versus baryonic mass m^, for the 
model A(ip) = exp(— 3ip 2 ). The solid line corresponds to the maximum value of <po allowed by so- 
lar-system experiments, and the dashed lines to (fo = ("zero- mode"). The dotted lines correspond 
to unstable configurations of the star. 




FIG. 2. Critical baryonic mass m cr versus the curvature parameter [3 within the quadratic 
models A{ip) = exp(i/3(/? 2 ). 
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FIG. 3. Dependence upon the baryonic mass rfiA of the coupling parameters a a, Pa, the inertia 
moment I a, and its derivative Slnl^/c^o- These plots correspond to the model A(ip) = exp(— 3ip 2 ) 
and the maximum value of cpo allowed by solar-system experiments. As in Fig. El the dotted lines 
correspond to unstable configurations of the star. 



Ka = -a A (3hiI A /a(po) 




FIG. 4. Parameter = —aA{dlnlA/dtpo) versus the Einstein inertial mass mAi within the 
model A(ip) = exp(— 3ip 2 ). The solid line corresponds to the maximum value of ipo allowed by 
solar-system experiments, and the dashed line to a ten-fold smaller value of (po (i-e, a 100 times 
smaller value of the Eddington parameter 7Edd — !)• 
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FIG. 5. The (^-7-^)1913+16 test for general relativity (GR), the Jordan-Fierz-Brans-Dicke 
theory (JFBD), and the quadratic model A(tp) = exp(+3<^ 2 ) [corresponding to a positive curvature 
parameter j3 = +6]. The widths of the three P lines correspond to l-a standard deviations. The 
(Jj th = (Jj cx p and 7 th = 7 exp lines are wider than l-a errors, and cannot be distinguished for the 
three theories. 




FIG. 6. The (w-7--P&)igi3+i6 test for the quadratic model A(ip) = exp(^/3(/9 2 ) with = —4.5, 
when ipo takes the maximum value allowed by solar-system experiments (left panel) , and a 30 times 
smaller value (right panel). In this Figure and the following ones, l-a deviations are smaller than 
the width of the lines. 
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FIG. 7. The (to-^-Pb) 1913+16 test for the quadratic model A(ip) = exp(— 3(/? 2 ) [i.e., (3 = —6], 
when ipo takes the maximum value allowed by solar-system experiments (left panel), and a 100 
times smaller value (right panel). 



m B /m Q 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 - 1 1 1 1 1 1 1 1 1 1 > m A /m 
0.5 1 1.5 2 2.5 

FIG. 8. The (07-7-3)1534+12 test for general relativity (GR), and for the quadratic model 
A{ip) = exp(— 3ip 2 ) [i.e., (3 = —6] when (po takes the maximum value allowed by solar-system 
experiments. The widths of the strips correspond to 1-a standard deviations. The arrow indicates 
the intersection of the three strips in general relativity. In the model (5 = —6, the three strips do 
not intersect. 
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FIG. 9. Regions of the (ao, (3o)-plane allowed by solar-system experiments and three bi- 
nary-pulsar experiments. In view of the reflection symmetry ao — > — «o, we plot only the upper 
half plane. The region allowed by solar-system tests is below the solid line. The PSR 0655+64 
test constrains the values of Qo and (3$ to be between the two dashed lines. The region allowed by 
the PSRs 1913+16 and 1534+12 tests lies to the right of the (approximate) wavy line. The region 
simultaneously allowed by all the tests is shaded. 
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TABLES 



TABLE I. Critical baryonic mass m cr (and critical Einstein mass m CT ) for some values of the 
curvature parameter (3 within the quadratic models A(ip) = exp(i/3( ) 5 2 ). 



p 


m cr /m 


m cr /m 


-10 


0.69 


0.66 


-9 


0.78 


0.74 


-8 


0.89 


0.84 


-7 


1.04 


0.98 


-6 


1.24 


1.16 


-5.5 


1.38 


1.28 


-5 


1.56 


1.43 


-4.5 


1.84 


1.65 


-4.35 


2.01 


1.78 
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